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ABSTRACT 

The location accuracy of the BeppoSAX Wide Field Cameras and acute 
ground-based followup have led to the detection of a decaying afterglow in 
X rays and optical light following the classical gamma-ray burst GRB 970228. 
The afterglow in X rays and optical light fades as a power law at all wave- 
lengths. This behaviour was predicted for a relativistic blast wave that radiates 
its energy when it decelerates by ploughing into the surrounding medium. Be- 
cause the afterglow has continued with unchanged behaviour for more than a 
month, its total energy must be of order 10 51 erg, placing it firmly at a red- 
shift of order 1. Further tests of the model are discussed, some of which can 
be done with available data, and implications for future observing strategies 
are pointed out. We discuss how the afterglow can provide a probe for the 
nature of the burst sources. 

Key words: gamma-rays: bursts - X-rays: transients - optical: transients - 
stars: neutron 



1 INTRODUCTION 

Gamma-ray bursts are such a mystery in large part because in the more than three decades 
since the first report on them (Klebesadel, Strong, & Olson 1973) they have remained in- 
visible in any radiation other than X and 7 rays. This situation ended in the early hours 
of February 28, when the Gamma- Ray Burst Monitor (40-1000 keV) and one Wide Field 
Camera (2-32 keV) on board the Italian-Dutch Satellite per Astronomia a Raggi X (Bep- 



© 0000 RAS 



2 Wijers, Rees, and Meszaros 

poSAX) triggered on a moderately bright gamma-ray burst and imaged it in X rays to give 
an error box only 6 arcmin across (Costa et al. 1997). It was still detected at a much fainter 
level by the MECS and LECS X-ray telescopes on BeppoSAX 8 hours after the burst, and 
again 3 days later, tightening the error box radius to under an arcmin. Optical images taken 
with the William Herschel and Isaac Newton Telescopes at La Palma starting only 20 hours 
after the burst reveal a fading source that is undoubtedly associated with GRB970228 (Van 
Paradijs et al. 1997), and therefore its location is now known with sub-arcsecond precision. 

Initial reports of an underlying host galaxy of roughly R = 24 (Van Paradijs et al. 1997) 
appear inconsistent with later HST images that show only a very faint, if any, underlying 
object (Sahu et al. 1997a,b). The issue of a host and its implications for the distance scale 
therefore remain to be resolved, and we shall be primarily concerned here with the evolution 
of the gamma-ray burst itself. We shall assume that its contribution to all optical detections 
thus far is not large. 

In Section ^| we briefly review the basic properties and predictions of a blast wave. In 
Section [| we discuss the light curves of GRB 970228 and GRB 970402 and their afterglow and 
show that they agree very well with the predictions of the fireball model (Rees & Meszaros 
1992, Meszaros & Rees 1997a). We then briefly discuss the distance scale in relation to our 
results and some possible complications due to less simple blast wave models (Sect. |j). We 
summarise our findings and discuss some implications for observing strategies in Section [|. 

2 THE SIMPLEST BLAST WAVE AND REMNANT EVOLUTION 

The simplest fireball remnant model is given by the diminishing emission of a forward blast 
wave moving ahead of a fireball, which continues to plough into an increasing amount of 
external matter beyond the deceleration radius at which the bulk Lorentz factor T first 
dropped substantially giving rise to the GRB (Rees & Meszaros 1992, Vietri 1997). Beyond 
this, the Lorentz factor decays as T oc r~ 3 / 2 oc t -3 / 8 , where t is the observer- frame time 
(Meszaros & Rees 1997a), and the spectrum is due to synchrotron radiation, from electrons 
accelerated to a power law N^) oc 7~ p above the minimum energy 7 m oc T imparted by the 
shock. The comoving intensity or energy spectrum is 




for v < u. 



for v > v x 



in 



(1) 
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where the synchrotron break frequency u m oc YB' r y^ a) and j3' = (1 —p)/2. Observationally 
(Band et al. 1993), a' = a+1 is around (single B' synchrotron would give 1/3 but sampling 
a range of B' and 7 m can easily give a range around that), while j3' — (5 + 1 is around — 1 
(ranging approximately from —0.5 to —2). This is the canonical GRB spectrum in 7-rays, 
the observed hu m being typically in the range of 50 keV to 2 MeV. 

The comoving-frame equipartition magnetic field is B' oc Y ~ t~ 3 / 8 , so the break fre- 
quency drops in time as is m oc TB'^ oc T 4 oc t~ 3 / 2 . At the same time, since the co moving 
electron density n' c oc Y and the comoving width of the emission shell AR ~ r/Y ~ t 5 / 8 , the 
comoving intensity V v oc n' e B' 2 ^AR/ '(B'^) oc n' c B'AR oc t -1 / 8 , so the observed flux as a 
function of observer time is (Meszaros & Rees 1997a) F Uin oc t 2 Y 5 I' Uin oc i° ~ constant. The ob- 
served break frequency z/ m crosses the X-ray or optical band at a time tx, op ~ (^/^x^p) 2 ^ 3 ^ 
after the break frequency lies at z/ 7 . For /iz/ 7) x, op = 100 keV, 5keV , and 2.3 eV, respectively, 
one has tx = 7t 7 and t op = 170t 7 . 

If the time-averaged flux at 7 rays of the GRB was F 7 , then for t < tx, op we have 
Fx,op — ^(^x.op/^m)"', After z/ m crosses the X-ray or the optical band, we have 

^x.op oc F Vm (u Xj0p /u m f oc t^ 13 ' = t 5 , for i > t x ,o P (2) 

The expansion Y oc r~ 3//2 oc t~ 3 / 8 lasts until the remnant becomes nonrelativistic, af- 
ter which the remnant enters a Sedov- Taylor phase. Here r oc t 2//5 , B' oc t~ 3 / 5 , 7 oc 
t" 6/5 , n' e AR oc t 2/5 and z/ m oc t~ 3 , F Uin oc t 3/5 , so that the optical (and X-ray, if it were 
detectable) flux has v > v m and goes as F op oc t( 3 + 15 / 3 ')/s ^ t" 12/5 for (3' ~ -1. Such a break 
in behaviour would occur when the blast wave has swept up a rest mass energy equal to its 
initial energy, at time 

= 1/3 ^ / M 1 t d / M 1/3 

y47rnm p c 5 y \ n -i / \ n -3 J 

where E is the initial explosion energy and n is the density of the surrounding medium, and 

the two numerical expressions are scaled to typical cosmological and halo cases, respectively. 

Such a break has, so far, not been seen in GRB 970228. Since it has now been followed for 

over a month, this firmly places it well beyond the halo of our Galaxy, if the blast wave 

model is indeed valid. 



3 CONFRONTATION WITH THE DATA 
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Table 1. Fluxes of GRB970228 and its fading counterpart. Time is measured from the burst trigger (February 28.124 UT). 
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(1) Costa et al. 1997a; (2) Palmer et al. 1997; (3) van Paradijs et al. 1997; (4) Matz et al. 1997; (5) Costa et al. 1997b; (6) 
Galama et al. 1997; (7) Margon et al. 1997; (8) Metzger et al. 1997b; (9) Smith et al. 1997; (10) Yoshida et al. 1997; (11) Groot 
et al. 1997b; (12) Groot et al. 1997c; (13) Tanvir & Bloom 1997; (14) Klose et al. 1997; (15) Sahu et al. 1997; (16) Sahu et al. 
1997; (17) Soifer et al. 1997; (18) Metzger et al. 1997a. 



3.1 Data on the light curve 

In Table |T] are listed all the fluxes from gamma-ray to radio that have been reported since the 
initial trigger. All upper limits are 3er. Optical and near-infrared magnitudes were translated 
into fluxes assuming they were all calibrated on the Vega system. Of course, since many 
reports are preliminary, the final calibrated values may differ somewhat from the table 
values, but since our emphasis is on trends of the fluxes when they change by a few orders 
of magnitude, these corrections will not affect the results reported here. The foreground 
reddening at the location of the burst is EiB — V) = 0.14 (Burstein & Heiles 1982), so the 
magnitudes were de-reddened using A B = 0.8, Ay = 0.4, A R = 0.3, Aj = 0.2, Aj^ H K = 0. 
The X- and 7-ray data need a somewhat more careful treatment, since they are taken 
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over much broader ranges of photon energy. TGRS reported clear evidence for a break in the 
spectrum, at E — 100 — 150 keV if an exponential cutoff was used (Palmer et al. 1997). On 
the other hand, the GRBM on SAX reported significant flux above 600 keV, so the spectral 
change is more likely to take the form of a power law spectrum that changes slope at some 
break energy E m , like the model of Band et al. (1993). The four fluxes reported for the 
initial strong 4-s spike (which appears contain the bulk of the burst fluence) are (in units 
of 10~ 6 ergcm -2 s -1 , with subscripts giving the energy range in keV): -F40-200 = 1-7 (Palmer 
et al. 1997), -Fi. 5 _7. 8 = 0.1, F 40 _6oo = 4, -F40-1000 = 6 (Van Paradijs et al. 1997). Assuming a 
spectrum that is flat below hu m = E m (a' = 0) and has slope (3 1 above it (eq. |l|), we can try 
to see which (/?', E m ) match the fluxes best. E m must be greater than about 40keV and j3' 
in the range —1 to —0.5 to get a satisfactory fit. Since we know from the TGRS spectrum 
that E m < 150 keV, this pins the parameters down reasonably well. To specify the X- and 
7-ray fluxes, we use E m = 40keV and f3' = —0.8. The gamma-ray flux is then given as the 
fitted value at 100 keV, -Fiookev = 4200 /iJy. The OSSE upper limits were likewise translated 
into values at 100 keV. 

The initial X-ray flux also follows from the spectral fit, and all other reported X-ray 
fluxes are translated into fiducial F u values at 5 keV as well. For X rays, this is unlikely to 
introduce an error of more than 30%, for gamma rays it may be a bit more. Accounting 
for uncertainties in reddening and the preliminary calibration, most optical points should 
have errors under 30%. This could affect the flux offset between light curves at different 
bands by up to that amount, but is unlikely to affect the inferred rate of decay in any single 
band by much, since the reddening corrections do not affect them. Moreover, in cases where 
more than one measurement in the same band by the same instrument is available, which 
presumably suffer from calibration effects in the same manner, the decay rate from that 
subset is quite consistent with that of all data in that band. The flux history of GRB 970228 
at various photon energies is plotted in Fig. [I]. 

3.2 Comparing GRB 970228 and 970402 with the blast wave model 

Fig. [I] clearly confirms the prediction that the flux of GRB 970228 should decline as a power 
law, in bands where more than 2 measurements well separated in time are available (X-ray, 
V, I). Moreover, a fit with free slope to data in those bands shows that the exponent of 
the decay, 5, is the same for all, as the model demands: 5 = —1.2. But then we predict 

© 0000 RAS, MNRAS 000, 000-000 



6 Wijers, Rees, and Meszdros 
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the slope of the spectrum, (3' = 25/3 = —0.8. The spectrum in the decaying part of the 
light curves (i.e. above the break) also follows from the flux ratios at fixed time, i.e. from 
the vertical offsets between the fitted lines. This independent measurement of the spectrum 
gives f3' = —0.78, remarkably close to the prediction. 

Since after fixing the vertical level of the fit to the V data all other levels are fixed in the 
model, we can show the predictions that follow from it for the other frequencies (Fig. [I]). The 
agreement is quite good all the way from X rays to K band. It is too early to say whether the 
occasional exception (e.g. the first J point) constitutes interesting extra behaviour. Overall, 
the agreement in such detail can be regarded as a strong confirmation of the model. 

It is undecided whether the initial gamma-ray burst derives from the same mechanism as 
the afterglow, i.e. is the initial part of the blast wave deceleration, or has a separate origin, 
e.g. internal shocks in the relativistic wind (Paczyhski & Xu 1994, Rees & Meszaros 1992) 
before the deceleration by the external medium begins. There is some indication in the 7-ray 
light curve of GRB 970228 that speaks against the former option. The leftmost dashed curve 
in Fig. [1] is the prediction for the 7-ray flux assuming it originates from the blast wave itself 
(and can therefore be predicted from the V and X-ray light curve). It only just falls below 
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the OSSE upper limits, and even then the break happens too early, as evidenced by the fact 
that the initial 7-ray flux lies above the curve. It is therefore possible that the gamma rays 
have a different origin. In that case, we can no longer rely on the initial 7-ray flux to set 
the level of emission from the blast wave before the break. This implies that the initial V 
magnitude of the burst could have been rather low, limited only by the first measurement 
(V = 20.9), rather than a (model dependent) value related to the initial 7-ray flux. 

The sparser data on GRB 970402, which was triggered on April 2.93 UT (Feroci et al. 
1997) also have some bearing on this matter. This burst lasted about 100 seconds with 
apparently no obviously dominating initial peak, and its peak flux was about 10 times lower 
than GRB 970228. It is the second GRB of which afterglow was detected, again 8 hours after 
the trigger, but a second attempt after 1.7 days only found an upper limit (Piro et al. 1997). 
Followup with 1-m telescopes at Siding Spring and SAAO (the burst is at declination —69°) 
0.7-1.9 days after trigger place upper limits of R = 21 on any fading counterpart (Groot 
et al. 1997a). From the X rays, and initial gamma rays, we once again infer that this burst 
has roughly the canonical spectrum of slope below and slope —1 to —0.8 above the break, 
so we expect it to fade in time as a power law with slope 5 = —1.5 to —1.2. If we assume 
all emission is from the blast wave, and hence the flux at all wavelengths starts at the same 
value as the initial 7-ray flux, then the time at which the break lies at 5keV follows from 
the X-ray detection: £5 = 4-40 s. The break should then arrive at R after tn =800-8000 s, 
consistent with the R limit. However, the predicted time for the gamma-ray break becomes 
^ioo =0.5-5 s, rather early for a burst that lasts over 100 s and is not dominated by an 
initial spike, as was GRB 970228. This discrepancy in the 7-ray timing, similar to the case of 
GRB 970228, may lend some support to the independent origin of the initial 7-ray emission. 

To distinguish between these cases, measurements in the first 15 minutes in X rays or in 
the first few hours in optical are needed. This could be achieved with present instruments 
because positions with 10 arcmin accuracy become available from the Wide Field Cameras 
on board BeppoSAX on a 1-3 hour time scale, and a BACODINE-type alert system could 
spread this knowledge to observatories around the world in seconds. Also, since re-pointing 
the satellite is not done until after some time, the WFC data immediately following the 
main burst could be a sensitive diagnostic of what the initial level of the X-ray afterglow 
might be. 
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4 DISCUSSION 

4.1 Galactic models 

Since they have much less energy than the cosmological ones, conventional matter-dominated 
fireball blast waves should fade in about a day (Sect. |2|), and thus cannot give rise to the 
observed long-term afterglow. 

In a galactic-halo model, the X-ray emission could be due to emission from a surface 
initially heated to T > 10 7 K, which might cool as a power law in time (regulated by some 
diffusion time dependent on the depth of deposition of the energy). The optical emission is 
much larger than the F op ~ 10~ 4 Fx expected from the Raleigh- Jeans tail of such a quasi- 
thermal surface emission. It is conceivable that the optical could be a separate component, 
perhaps associated with reprocessing of harder radiation by circumstellar material, in which 
case one might expect the optical first to increase, peak, and then fade, as the irradiating 
"hard" photons first cool to become UV and then become themselves optical, dropping off 
exponentially after that. 

The Comptonisation model for gamma-ray bursts in the Galactic halo (Liang et al. 1997) 
predicted the optical brightness 20 hr after the burst to be V = 25 from the observed X-ray 
flux (Smith et al. 1997), 4 magnitudes fainter than observed. Hence it is inconsistent with the 
behaviour of GRB 970228. This is due mainly to the fact that it predicts that the spectrum 
below the break fades rapidly as well, in contrast with the blast wave model. 

The radiation efficiency in galactic halo blast wave models would generally be low, 10~ 3 
or less, due to the smaller magnetic field strengths at a given time and the shorter expansion 
time scales at a given T. Thus one might expect faster turnoff of both the X-ray and optical 
emission, due to a decrease of the optical radiation efficiency. This would be in addition to 
the steepening caused by moving into the Sedov- Taylor nonrelativistic phase. The late-time 
radio afterglow should come from a large enough region to be resolved by VLBI. 

For sources within a few hundred kpc, one could expect the radio emission to turn on 
much sooner than in cosmological models, since it is quicker in reaching a given angular 
diameter, when the remnant still has sufficient internal energy to produce observable emis- 
sion. Furthermore, a galactic halo neutron star would not be disrupted (since it needs to 
repeat many times) and thus could, in principle, produce also coherent radio emission by 
pulsar-type mechanisms in the magnetosphere. 
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4.2 Less simple fireball and remnant evolution 

The decay model discussed in Sect. |2] is for the simplest external blast wave scenario, valid 
for a spherically expanding fireball or a flow which is uniform inside channels of smaller 
solid angle. Other types of emission have been also considered, e.g. from the reverse shock 
moving into the ejecta, or from ejecta with a frozen-in magnetic field, or from internal shocks 
(Meszaros & Rees 1997a). These decay more rapidly than F oc t( 3 / 2 )^' (eq. |2|). However, there 
are additional considerations which can lead to different decay laws, including some which 
are slower than t( 3 / 2 )^'. It is, for instance, probable that the high-I" fireball that gives rise to 
the GRB may be beamed (Meszaros & Rees 1992, 1997b), and in this case it is natural for 
T to be a function of angle, tapering off towards the edges of the beam. In a cosmological 
compact binary disruption scenario, a central object (presumably a black hole) forms, which 
is surrounded by a temporary torus of neutron star debris. A very high bulk Lorentz factor 
jet emerges around the rotation axis, which would taper off into a slower wind, as it mixes 
with the increasingly baryon loaded wind at large angles which must result from the super- 
Eddington radiation from the torus. As a simple example, we can take a dependence T oc 9~ k 
(Meszaros & Rees 1997b), where 9 is the angle from the axis of rotation, outside of the main 
part of the jet. The deceleration of the fireball commences at a radius oc r _8//3 , which is 
larger for larger 9. Since radiation can be seen by the observer from within a cone of angle 
~ r _1 , the radiation from increasingly larger angles will be detectable by the observer, at 
later times which are given by the deceleration time t&{8) oc r~ 8 / 3 for that Y ~ 9~ k . This 
time is of order a day for V ~ 10, and a month for r ~ 3. In the simplest case where 
the burst progenitor supplies equal amounts of energy into equal logarithmic intervals of 
9, the bolometric flux detected as a function of time will then be F(t) oc E(9)/t oc £ -1 . 
If the spectrum is again of the form F u oc [ u a \ u 13 ' ] below and above a time dependent 
break v m , and we take the field to be a constant fraction of the equipartition value, then 
u m oc YB'^ 2 oc T 4 oc t~ 3 / 2 , and F Vm ~ L/u m oc t~H 3 ' 2 oc t 1//2 . If the radiative efficiency does 
not change in time and the spectrum evolves homologously, then we have 
' Fvm {!^Y cctWocf 1 / 2 , 

for t < t x ,o P , if a' ~ 

Fx op ~ s (4) 

for t > i x ,op, if P' ~ -1 
If the optical radiative efficiency varies in time, the simple homologous behaviour above 
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introduces additional changes. If we again assume that fields are a constant fraction of 
equipartition, then v oc TB'^ 2 oc T 2 7 2 , and the electrons responsible for optical radiation 
must satisfy -y op oc T -1 , so the ratio of comoving expansion and cooling times is t' cxp /t' cool oc 
T~ 2 / 3 oc t 1 / 4 . Then the above time power laws remain the same as long as the optical efficiency 
£ op = max[l, (texpAcooi)] i s unity, or they could be flatter by t 1//4 if the optical efficiency was 
initially lower than unity and grows (e.g. the higher energy electrons emitted more efficiently 
at higher photon energies than at optical). There are, of course, other possibilities; for 
instance, if the field is frozen-in, rather than turbulent, a different B 1 dependence needs to 
be used, etc. The above argument, however, indicates that there are ways in principle to 
explain even slower decay time scales than the simple oc £( 3 / 2 )^', as well as faster ones. The 
afterglow might then probe the geometry of the emission, and for suitably oblique viewing 
angles one would see afterglow without a burst. 

It is important to note that these are testable differences, because the dependence of 
the power-law slope of the temporal evolution, 5, on the slope of the spectrum, /?', is not 
the same for eqs. |2] and |j. In the case of GRB 970228, with f3' = —0.8, the predictions are 
5 = —1.2 for the simple case and 5 = —0.7 from eq. |], so the simple case is favoured for it. 

5 CONCLUSION 

As expected, the first detection of a gamma-ray burst in the optical has greatly furthered our 
understanding of these enigmatic objects. We have found that the simplest fireball model for 
a gamma-ray burst and its afterglow agree very well with the data obtained for GRB 970228 
and GRB 970402. The longevity of the afterglow argues strongly, within the context of 
that model, that GRB 970228 occurred at substantial cosmological redshift. This raises the 
interest in unveiling the nature of the faint extended object coincident with the fading burst. 
Its faintness would suggest that if it is a galaxy it probably has z > 1, whereas the burst itself 
was moderately bright. For a no-evolution, standard candle cosmological burst distribution, 
its redshift as inferred from the peak flux in 7 rays would be smaller, perhaps z ~ 0.3. This 
would argue either for a modest width of the GRB luminosity function, or for significant 
evolution of their rate density. Neither is implausible, e.g. the star formation rate evolves 
so strongly with redshift (Lilly et al. 1996) that if the GRB rate were proportional to it 
(which a merging neutron star scenario would naturally demand) the GRB log N — log P 
distribution would seem Euclidean to z = 1. 
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We note that optical flashes perhaps even brighter than V ~ 19 — 20 from GRB at 
cosmological distances (Meszaros & Rees 1997a) would be of enormous importance as a 
tool to study absorption lines from the early intergalactic medium, features due to the ISM 
in the host galaxy etc. (Miralda-Escude, Rees & Meszaros 1997), particularly if they arise 
from redshifts as large as z ~ 3 — 5. Even at more moderate redshifts, a good spectrum 
obtained with a large telescope in the first day could settle distance scale definitively by 
showing redshifted absorption lines or even a Lyman alpha forest^; even better would be a 
space-based instrument with access to the UV part of the spectrum. 

We also argue that further significant tests of afterglow models depend crucially on 
catching the fading counterpart within the first 0.5 to 3 hours after the trigger. This is 
feasible in principle with BeppoSAX since the first (10 arcmin precision) position from the 
WFC can be derived within the hour. The history of the first two afterglow detections show 
that important information comes from upper limits as well, especially at early times. So 
even with modest telescopes rapidly following up a GRB alert is very worth while. Since the 
light curve appears to evolve as a power law, one should try to integrate no longer than the 
time that passed since the GRB trigger in very early measurements, since the brightness of 
the transient may change on that time scale. 
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